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1 Introduction 

Let W be a closed subspace of a separable complex Hilbert space and let I be a finite or countable 
infinite set. A frame T = for W allows for linear (typically redundant) and stable encoding¬ 

decoding schemes of vectors (signals) in W (see Section [2.11 for definitions and technical results). 
Indeed, if V is a closed subspace of Ti such that V 0 W"*" = % (e.g. V = W) then it is possible to 
find frames Q = {fl’diei for V such that 

/ = fi , for / E W. 

iei 

The representation above lies within the theory of oblique duality (see [T71 [181 [22l [23]). fo applied 
situations, it is usually desired to develop encoding-decoding schemes as above, with some additional 
features. In some cases, we search for schemes with prescribed properties (e.g., for which the 
sequence of norms {||/i|p}iGii as well as the spectral properties of the family E are given in advance) 
leading to what is known in the literature as frame design problem (see 13 0 [mini [351 [12]). In 
other cases, we search for numerically robust oblique dual pairs {T, Q) as above, leading to what is 
known as optimal frame designs (13 m [23 [33 US HZ]). 

In their seminal work [3, Benedetto and Fickus introduced a functional defined on finite sequences 
of (unit norm) vectors E = (where = {1,... , n}), the so-called frame potential, given by 

FP(.F)= . (1) 

j Gin 

In case dim 77 = d E N then one of their major results shows that tight unit norm frames can 
be characterized as (local) minimizers of this functional, among unit norm frames. Since then, 
there has been interest in (local) minimizers of the frame potential within certain classes of frames, 
since such minimizers can be considered as natural substitutes of tight frames (see for example 
mm)- Recently, there has been interest in the structure of minimizers of other potentials such 
as the so-called mean squared error (MSE) (see [2411371 H2] and the references therein). Both the 
frame potential and the MSE are examples of the so-called convex potentials introduced in [33 • It 
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turns out that minimizers of these convex potentials share the spectral and geometrical structure 
of minimizers of the frame potential. Now, it is a well known fact that in case V = W then tight 
frames J- for W - i.e. minimizers of convex potentials - give rise to optimal (numerically robust) 
dual pairs (J^, Q\ Therefore, it seems apparent that in the general case VSkV"*" the construction of 
robust oblique dual pairs (J-”, Q) is related with the construction of frames T which are minimizers 
of convex potentials (e.g. the frame potential). 

It turns out that there is a single notion that ties all the previous problems together namely, the 
majorization preorder. Indeed, majorization is the key notion behind the frame design problems 
(see [anaiii]) through natural extensions of the Schur-Horn theorem from matrix analysis (see 
[3 El dm EH). Moreover, the relation between majorization and tracial inequalities with respect to 
convex functions allows to apply this notion in the study of convex potentials ([3 ESI EH ESI SO]). 
Unfortunately, the convex potentials considered in [36] (in particular, the frame potential) can only 
be defined for finite frames. Hence, in the infinite dimensional context we loose a tool which have 
proved useful as a measure of stability for frames in finite dimensional Hilbert spaces. 

In this paper we show that there are natural analogues of the convex potentials (and in particular, 
of the frame potential) in the context of Bessel sequences of integer translates E{J^) of a finite 
family of vectors T = that lie in a finitely generated shift invariant (FSI) subspace W of 

L^(M^). We take advantage of the detailed structure of E{T) as a fibered family over the torus 
(see [71 El US) and the theory of range functions for shift invariant (SI) subspaces and introduce 
the potential {E{E)) associated to a convex function cp and W as an integral over of the 
corresponding potentials on the fibers (for related approaches to different problems in SI subspaces 
see [D ESI EH]). In order to verify that our definition is a natural extension of the convex potentials 
for finite frames we show that under natural normalization conditions, a family of integer translates 
E{P) that is a tight frame for a FSI subspace W is a minimizer of the convex potential associated 
to if and W. 

The convex potentials in FSI subspaces raise several questions related with optimal design problems. 
In particular, given FSI subspaces W, V such that VQW"*" = L^(]R^) and a finite family E = 
such that E{E) is a frame for W, we consider the problem of designing optimal oblique duals E{Q) 
which are translates of a family Q = in V and such that Q satisfies the norm restrictions 

WdiW'^ — w > 0. In order to deal with this problem we develop two new tools in the 

context of frames of translates. On the one hand, we obtain what we call the fine spectral structure 
of shift generated (SG) oblique dual frames of the fixed frame E{E), which is a detailed description 
of the eigenvalues of the measurable field of positive operators defined on corresponding to the 
frame operators of SG oblique duals of E{E). As a consequence, we derive necessary and sufficient 
conditions for the existence of a tight SG oblique dual E{Q) of E{E). On the other hand, we consider 
the water-filling construction (both for functions in probability spaces as well as for measurable field 
of positive finite-rank operators) and show that this construction leads to optimal solutions of the 
oblique dual design problem; this is achieved by showing that the water-filling constructions are 
optimal with respect to majorization (considered in the general context of probability spaces) which 
is a result of independent interest. With these tools we completely solve the problem of designing 
optimal oblique dual frames with norm restriction mentioned before; it turns out that these optimal 
SG oblique duals are more stable than the so-called canonical oblique dual. We point out that the 
structure of the optimal solution is obtained in terms of a global analysis. As a byproduct we extend 
the so-called Fan-Pall interlacing theorem from matrix analysis to the context of measurable helds 
of positive matrices. 

The paper is organized as follows. In section [2] we describe several preliminary notions and facts 
from frame theory, SI subspaces, frames of translates and majorization theory in probability spaces. 
In section El we introduce the convex potentials for frames of translates and show that are natural 
extensions of the convex potentials for finite frames. In section E] we recall several facts from the 
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theory of oblique duality in FSI subspaces and obtain the precise value of the aliasing norm corre¬ 
sponding to the consistent sampling in this setting. Then, we describe the hne spectral structure of 
oblique duals of a fixed SG frame. Since this result depends on an extension of the Fan-Pall interlac¬ 
ing theory, its proof is presented in an appendix (see section E]). In section 0 we study the problem 
of optimal design of oblique dual frames E{Q) - of a fixed finitely SG frame E{E) - which satisfy 
certain norm restrictions. We first show that the water-filling construction for positive functions 
in probability spaces is optimal with respect to sub-majorization within a natural set of functions. 
We then construct optimal SG oblique duals with norm restrictions and explain the relation of our 
construction with a natural (non-commutative) water-filling construction for measurable field of 
positive finite-rank operators. The paper ends with an appendix section in which we develop the 
Fan-Pall interlacing theorem for measurable fields of positive matrices as well as some consequences 
of this result. 

2 Preliminaries 

In this section we recall some basic facts related with frame theory, oblique duality and shift 
invariant (SI) subspaces of L^(M^). At the end of this section we describe majorization between 
functions in arbitrary probability spaces. 

General Notation 

Throughout this work we shall use the following notation: the space of complex d x d matrices 
is denoted by Alrf(C) and A4d(C)“'" denotes the set of positive semidefinite matrices. Ql{d) is the 
group of invertible elements of and Ql{d)~^ = A4d(C)“'“ fl Ql{d). If T € A4^{C), we denote 

by ||r|| its spectral norm, by rkT = dimi?(T) the rank of T, and by trT the trace of T. 

Given d E N we denote by = {1,..., d} C N. For a vector x E we denote by (resp. x^) the 
rearrangement of x in decreasing (resp. increasing) order. We denote by (M'^)'*' = {a: E : x = x^'\ 
the set of downwards ordered vectors, and (M'^)^ = {x E : x = x^}. 

Given S E AIrf(C)+, we write A(S) = X^{S) = (Ai(S'), ... , Arf(S')) E (M'^)'*' the vector of eigenvalues 
of S' - counting multiplicities - arranged in decreasing order. Similarly we denote by A^(S') E (M'^)^ 
the reverse ordered vector of eigenvalues of S. 

If VF C is a subspace we denote by Pw G the orthogonal projection onto W. Given 

X, y E we denote by x ( 8 > y E Alrf(C) the rank one matrix given by 

X ® y {z) = {z , y) X for every z E . (2) 

Note that, if x 7 ^ 0, then the projection = T’spania:} = x ® x . 

2.1 Frames for subspaces and oblique duality 

In what follows Ji denotes a separable complex Hilbert space and I denotes a finite or countable 
infinite set. Let W be a closed subspace of Ti: recall that a sequence P = in W is a frame 

for W if there exist positive constants 0 < a < 6 such that 

a\\ff <b\\ff for every /eW. (3) 

iei 

In general, if P satisfies the inequality to the right in Eq. ([3]) we say that is a 6 -Bessel sequence. 

Given a Bessel sequence P = {fi}i£i we consider its synthesis operator Tjr e given by 

lV((«i)iGi) = Yli&P fi which, by hypothesis on P, is a bounded linear transformation. We also 
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consider G called the analysis operator of F, given by 7>(/) = and the 

frame operator of F defined by Sjr = TjrTfr- It is straightforward to check that 

{S^f, f) = '^\ if, fi) 1 ^ for every f eV. . 

iei 

Hence, Sjr is a positive semidefinite bounded operator; moreover, a Bessel sequence in W is a 
frame for W if and only if Sjr is an invertible operator when restricted to W or equivalently, if the 
range of Tjr coincides with W. 

In order to describe oblique duality, we fix two closed subspaces V, W such that W-*- © V = Ti, 
that is such that W"*" + V = T-L and W"*" fl V = {0}. Hence, W'^ is a common (algebraic) complement 
of W and V. It is well known that in this case -Pyvlv : V ^ W is a linear bounded isomorphism so, 
in particular, we see that dim V = dim W as Hilbert spaces. Moreover, the conditions W'^ © V = 
and W © V"*- = H are actually equivalent. 

Fix a frame F = {fi}i£i for W. Following [221E3] (see also [IS]), given a Bessel sequence Q = {gi}i£i 
in V we say that ^ is a (oblique) V-dual of F if 

f = 9i) fi for every / G W . 

i£l 

It turns out (see [22l[23|) that if ^ is a V-dual of F then Tjr Tg = Py^j , where Pyy /denotes 
the oblique projection with range W and null space V"*". In this case, by taking adjoints in the 
previous identity we also get that Tg Tfr = = -fy// w-l • Hence, Tg is onto V and then Q is 

a frame for V; moreover, we obtain the reconstruction formula 

9 = '^{ 9 , fi) 9i for every 5 G V . 
iei 

We consider the set of oblique V-duals of F given by 

Vy{F) = {Q = {9i}i& is a V-dual of 7'} . (4) 

Remark 2.1. Let F = {/ijigi be a frame for W. If we set V = W then a Bessel sequence ^ in W 
is a W-dual of F if it is a dual frame for W in W in the classical sense (see [TB|) i.e. Tg Tfr = Pw- 
Hence 

^w(-^) = ^(•^) '= {G = {9i}iei is a dual frame for 7" in W } . 

Recall that there is a distinguished (classical) dual, called the canonical dual of F, denoted F"^ = 
if*} jgi, given by ff = for i G I, where denotes the Moore-Penrose pseudo-inverse of the 
(closed range positive semidefinite operator ) Sjr. 

In the general context of oblique duality there also exists a distinguished V-dual for F, the so-called 
eanonical V-dual, which we denote by 

^V={f*i}i& given by f*.=Py//,^^f* for every iGl, 

where 7"^ = {ff'}i^i is the canonical dual of F. It turns out that the encoding-decoding scheme 
based on the oblique dual pair (7”, Fy ) has several optimality properties (see [221123] ) ■ A 
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2.2 Shift-invariant subspaces and frames of translates 

In what follows we consider (with respect to Lebesgue measure) as a separable and complex 

Hilbert space. Recall that a closed subspace V C L^(M^) is shift-invariant (SI) if / € V implies 
Tef G V for any £ G where Tyf{x) = f{x — y) is the translation by y € For example, if we 
take a subset A C L^(M^) then, 

<S(^) = spin {Tef :feA, (5) 

is a shift-invariant subspace called the SI subspace generated by A. Indeed, 5(>I) is the smallest SI 
subspace that contains A. We say that a SI subspace V is finitely generated (FSI) if there exists 
a hnite set A C L^(M^) such that V = S{A). In this case, the length of V is the smallest cardinal 
ff{A) such that S{A) = V. 

In order to describe the fine structure of a SI subspace we consider the following representation of 
L^(]R^) (see [3 El [271 US] and m for extensions of these notions to the more general context of 
actions of locally compact abelian groups). Let T = [—1/2,1/2) be endowed with the Lebesgue 
measure and let L^(T^,.^^(Z^)) be the Hilbert space of square integrable I’^(Z^)-valued functions 
that consists of all vector valued measurable functions (/ : —>■ t'^(Z^) with the norm 

Uf = \\Hx)\\%^zic) dx < oo. 

Then, F : f {ifi)) defined for / € L^R^) n L‘^{R^) by 

, rf{x) = if{x + £)fi^^,, (6) 

extends uniquely to an isometric isomorphism between L^(M^) and L^(T^, ^^(Z^)); here / denotes 
the Fourier transform of / G L^(M^). 

Let V C L^(M^) be a SI subspace. Then, there exists a function Jy : ^ { closed subspaces of 

^^(Z*^)} such that: if Pjy(x) denotes the orthogonal projection onto Jv{x) for x G T*', then for every 
y G £‘^{lfi) the function x i-)- {Pj^(x) ^ is measurable and 

V = {/ G L2(Mfc) ; r/(x) G Jv(x) for a.e. x G T^}. (7) 

The funcion Jy is the so-called measurable range function associated with V. By [51 Prop.1.5], Eq. 
([7|l establishes a bijection between SI subspaces of L^(]R^) and measurable range functions. In case 
V = 5* (.4.) C L2(R^) is the SI subspace generated by .4. = {/i* : i G 1} C L^(M^), where I is a finite 
or countable infinite set, then for a.e. x G we have that 

Jv(x) = {Fhj(x) : i G . ( 8 ) 

Recall that a bounded linear transformation S G L{Lfi{R^)) is shift preserving (SP) if 
for every ^ G Z^. In this case (see [HI Thm 4.5]) there exists a (weakly) measurable field of operators 
[5](.) : > £'^{lfi) (i.e. such that for every y G £'^{lfi) the function 9 x i-)- ([S']^;^, fi) is 

measurable) and essentially bounded (i.e. the function 3 x i->- || [5]^: || is essentially bounded) 
such that 

[5],(F/(x)) = F(5/)(x) for a.e. x G T^ / G L2(]Rfc) . (g) 

Moreover, ||5|| = esssup^.g-j-fc || [ 5 ] 2 ; ||. Conversely, if s : —)■ L{£‘^{7fi)) is a weakly measurable and 
essentially bounded field of operators then, there exists a unique bounded operator S G L{L‘^{R^)) 
that is SP and such that [S] = s. For example, let V be a SI subspace and consider P\> G L(L^(M^)), 
the orthogonal projection onto V; then, P\> is SP so that [Py] : —)• L{£‘^{7fi)) is given by 

[Pv\x = Pjv(x) i- 6 -) the orthogonal projection onto Jy(x), for a.e. x G T^. 
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The previous notions associated with SI subspaces and SP operators allow to develop a detailed 
study of frames of translates. Indeed, let T = {fi}i£i be a (possibly finite) sequence in 
We denote by E{T) the family of translates of T", namely E{E) = For a: G T^, let 

rT'(x) = {Tfi{x)}i^i which is a (possibly finite) sequence in £^(Z^). Then [8ll33] E{E) is a 6-Bessel 
sequence if and only if rj-'(x) is a 6-Bessel sequence for a.e. x € T^. In this case, we consider 
'^TT{x) ■ ^^(I) and 5'r:F(x) ^ —>■ £^(Z^) the synthesis and frame operators of rj^(x), 

respectively, for x € it is straightforward to check that Se{j^) is a SP operator. 

If T” = and Q = {gi]i^i are such that E{E) and E{Q) are Bessel sequences then (see [25103]) 

the following fundamental relation holds: 

\^E{g) ~ '^rg{x) '^vt{x) ’ x . ( 10 ) 

Eq. ([To]) has several consequences; indeed, if W is a SI subspace of L^(]R^) and we assume further 
that E, Q ^ then: 

1. For every f,g€^ L^(M^), 

{Se{E) f,g)= / (‘S'rj-(x) r/(x), Tg{x))p(^^k) dx . 

This last fact implies that = Syt{x) for a.e. x € T^; moreover, it also holds that 

E{E) is a frame for W with frame bounds 0 < a < 6 if and only if TE{x) is a frame for Jy\;{x) 
with frame bounds 0 < a < b for a.e. x € (see m)- 

2. Since [Pvv]x = Pj^ix) for ri-e. x G then E{Q) is a (classical) dual for E{F) in W if and 
only if F^(x) is a (classical) dual for Fjr(x) in Jy\){x) for a.e. x G (see [51 l28l [29] i. 


2.3 Majorization in probability spaces 

Majorization between vectors (see mm) has played a key role in frame theory. On the one 
hand, majorization allows to characterize the existence of frames with prescribed properties (see 
o da E]). On the other hand, majorization is a preorder relation that implies a family of tra- 
cial inequalities; this last fact can be used to explain the structure of minimizers of the so-called 
Benedetto-Fickus frame potential dans]) as well as more general convex potentials for finite frames 
(see [35l [36l EH [Ml HQ])- ^ section we extend the notion of convex potentials to the con¬ 

text of Bessel families of translates of finite sequences; therefore, we will need the following general 
notion of majorization between functions in probability spaces. 

Throughout this section the triple (X, fb, g) denotes a probability space i.e. fb is a a-algebra of 
sets in X and g is a probability measure defined on X. We shall denote by = {/ G 

L°°{X,fi) : / > 0}. For / G L°°{X,fi)~^, the decreasing rearrangement of / (see |S]), denoted 
f* : [ 0 , 1 ) —>■ M+, is given by 

f*{s) = sup {t G : g.{x G X : /(x) > t} > s} for every s G [0,1). (11) 


Remark 2.2. We mention some elementary facts related with the decreasing rearrangement of 
functions that we shall need in the sequel. Let / G L°°{X,g)~^, then: 


1 . 

2 . 


f* is a right-continuous and non-increasing function. 


/ and /* are equimeasurable i.e. for every Borel set ^4 C M then g,{f~^[A)) = 
where |i?| denotes the Lebesgue measure of the Lebesgue measurable set R C M. 
implies that for every continuous (p : then: cpo f £ L°°{X, g) iff (/? o f* 

and in this case 


/ f dg= [ ipo f 

Jx Jo 


dx. 


In turn, this 

gl-([0,i]) 
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3. If 5 € L°°{X, fi) is such that f < g then 0 < f* < g*] moreover, in case f* = g* then f = g. 

4. If we consider the probability space ([0,1], dt) - Lebesgue measurable sets in [0,1] with 
Lebesgue measure - then f* G L°°{[0, 1], dt) is such that (/*)* = f*■ 

5. If c G M is such that / + c > 0 then (/ + c)* = /* + c. A 

Definition 2.3. Let f,g (z and let /*, g* denote their decreasing rearrangements. We 

say that / submajorizes g (in {X,X^g)), denoted g -<w f, if 

ns ps 

/ g*(t) dt< f*{t) dt for every 0 < s < 1 . 

Jo Jo 

If we further have that g*{t) dt = f*{t) dt we say that / majorizes g and write g < f- A 

In order to check majorization between functions in probability spaces, we can consider the so-called 
doubly stochastic maps. Recall that a linear operator D acting on L°^{X,fj,) is a doubly-stochastic 
map if D is unital, positive and trace preserving i.e. 

D{lx) = lx, D(L°°(X,^)+) CL-(X,/r)+ and [ D{f){x) d^i{x) = [ /(x) d/r(x) 

Jx Jx 

for every / G L°°{X,fi). It is worth pointing out that D is necessarily a contractive map. 

Our interest in majorization relies in its relation with integral inequalities in terms of convex 
functions. The following result summarizes this relation as well as the role of the doubly stochastic 
maps (see for example mm)- 

Theorem 2.4. Let f,g& L°°{X,fi)^. Then the following conditions are equivalent: 

1 - 5 ^ /; 

2. There is a doubly stochastic map D acting on L°°{X,pl) such that D{f) = ( 7 ; 

3. For every convex function gp : ^ we have that 

ip{g{x)) dg{x) < [ g^{f{x)) dg{x) . (12) 


L' 


lx 


In case we only have g -<w f then Eq. m holds if we assume further that (p is an increasing convex 
function. □ 

Example 2.5. The operator D given by D{f) = (J”^ / dp) • lx is a doubly stochastic map. Hence, 
we get the majorization relation (J)^ / dp) • lx R /. Therefore, if p : M'*' —>■ M'*' is any convex 
function and / G L°°{X,p)~^ then, by Theorem 12.41 we have that 

[ f dp)= [ pH [ f dp) - lx{x)) dp{x) < [ pif{x)) dp{x) , ( 13 ) 

Jx Jx Jx Jx 

which is an instance of the classical Jensen’s inequality. Using the previous facts, notice that if 
c G M is such that 0 < c < f dp then it is easy to see that c ■ lx -<w f ■ A 

The following result will play a key role in the study of the structure of minimizers of -<w within 
(appropriate) sets of functions. 

Proposition 2.6 ([IS])- Let f,g& L°°{X,p)~^ such that g -<w f- If there exists a non-decreasing 
and strictly convex function p : R'*' —>■ R^ such that 


/ Ji’ifix)) dp{x) = [ p{g{x)) dp{x) then g* = f* 
Jx Jx 


□ 


With the notations of Example 12.51 notice that Proposition 12.61 implies (the known fact) that if p 
is strictly convex and such that equality holds in Jensen’s inequality Eq. (fl^ then f* = /x / dh 
and hence f = Jx f dp. 


3 Convex potentials for sequences of translates in FSI spaces 

We begin by describing the convex potentials for finite sequences of vectors with respect to a finite 
dimensional subspace. We consider the sets 

Conv(M"'') = {(p : —)• M'*' , (^ is a convex function } 

and ConVs(K'’') = {(/? G Conv(R+) , (p is strictly convex }. 

Now, given p G Conv(R+) and a finite dimensional subspace W C H, then the convex potential 
associated to {p, W), denoted by , is defined as follows: for a finite sequence P = G 

with frame operator Sjr g 


P^{P) = tT[piS^)Pw] ( 14 ) 

where p{Sjr) G is obtained by functional calculus and tr(-) denotes the usual (semi-finite) 

trace in L{'H). Notice that by construction, Py^ Sjr = Sj^ P\y = Sj^: then, it is clear that 

P^{P) = Y.p{h{Sr)). ( 15 ) 

i&d 


where d = dimW and (Ai( 5 jr))jgj^ g (R"'')'^ denotes the vector of eigenvalues of the positive 
operator G L(>V)'’“, counting multiplicities and arranged in non-increasing order (we use the 

convention Iq = 0)- In particular, if y? G Conv(R"’') is such that :^(0) = 0 we get that 

P^^(.F)=tr(:^(5^))=tr((^(Gj-)), 

where the n x n matrix Gjr = {{fi, fj))i,j£i„ is the Gramian matrix of the finite sequence P. That 
is, if (/?(0) = 0, then P^ = P^ does not depend on W. For example, in case p{x) = then 
P^(P) = Pip{P) coincides with the frame potential: indeed, by Eq. ([1]) we have that 

P^{P) = P^{P) = tr(53.) = tT{G%) = \{n, /,)| 2 = FP (.F) . (16) 

*> i sin 

For p G Conv(R'*') and a finite dimensional subspace W CP, P^(P) is a measure of the spread of 
the eigenvalues of the frame operator of P = {fi}i£i„ G That is, (under suitable normalization 
hypothesis on P) the smaller the value P^(P) is, the more concentrated the eigenvalues of 5^|w G 
are. This is the main motivation for considering these convex potentials (see [3611371 Hnll42j i. 

Next we extend the notion of convex potential to the context of finitely generated shift invariant 
systems in FSI subspaces. 

Definition 3.1. Let W be a FSI subspace in L^(R*'), let P = G be such that E{P) is 

a Bessel sequence and consider p G Conv(R+). Then the convex potential associated to {p,W) on 
E{P), denoted P^{E{P)), is given by 

P^{E{P)) = [ P;^^^^\TP{x)) dx (17) 


where PI^'^^^\tP{x)) = tr((/?(5rj-(3;)) [P'w]x) is the convex potential associated with {p,Jy\;{x)) of 
the sequence rT'(x) = {T fi{x)}i^i^ in f’^(Z^), for every x G A 

Next we develop some notions and tools in order to show that the right hand side in Eq. (HZD is 
well defined, namely that the function 9 x i—)■ PI^'^^^\tP{x)) is integrable. 
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Let T = {fi}m„ be a finite sequence in L^(M^) such that E{F) is a Bessel sequence. Recall that 
in this case Se(j^) is a SP operator and that for a.e. x G T^, = 5'rj-(a;) G L(f^(Z^))+ is a 

positive and finite rank operator. 

The next lemma is a reformulation of a result in |43j concerning the existence of measurable functions 
of eigenvalues and eigenvectors of measurable fields of positive semidefinite n x n matrices. 

Lemma 3.2. Let W be a FSI subspace in L^(R^) and let T = he such that E{F) is 

a Bessel sequence. Then, there exist: 


1 . 


a measurable function r : —>■ N>o and measurable vector fields Vj : —>■ for j G In 

such that r{x) < n and {vj{x)Y-^^\ is an orthonormal system in Jw{x) for a.e. x G T^; 


2 . 


bounded measurable functions Xj : — >■ R+ for j G In, such that Ai > ... > An , Aj(x) = 0 if 

j > r{x) and 

r(x) 

[5^(_7r)]a; = ^ Aj(x) Vj{x) (g> Vj{x) , for a.e. x G . (18) 


If we assume further that E{E) is a frame for W then r{x) = dim Jvv(a;) ond {vj{x)Y^Y\ is an 
orthonormal basis (ONB) for Jy\>{x) for a.e. x G 

Proof. Consider the measurable field of positive semidefinite matrices G : ^ A4n(C)^ given 

by the Gramian G{x) = {{Tfi{x),rfj{x)))ij^i^, for x G T^. Notice that G{x) is the matrix 
representation of G L(C"’)+ with respect to the canonical basis of C” for x G In 

particular, if b denotes a Bessel (upper) bound of E{E) then 


ess sup,^gT'=l|G'(a=)ll = ess sup^^gp/^ II = I|5’e(j^)II < b , 

by the remarks at the end of SectionWe set r(x) = rk(G(x)) = rk(S'r_p(j,)) for x G T^; therefore 
r(-) : ^ N>o is a measurable function such that r{x) < n for x G T^. Hence, by considering 

a convenient finite partition of into measurable sets we can assume, without loss of generality, 
that r{x) = r G N for a.e. x G . 

Using results from [43] . we see that there exist measurable functions Xj : —>■ R'*' and measurable 

vector fields Uj : ^ C", for j G In, such that: Aj(x) > Aj+i(x) for j G In-i, {'Uj{x)}j^i„ is an 

ONB of C"’ and G{x)uj{x) = Xj{x)uj{x) for j G In and a.e. x G T^. In particular, the functions 
Xj : ^ R+ satisfy 0 < Xj{x) < ||G(x)|| < h for a.e. x G T^, j G In; these remarks prove item [2] 

above. 

Take the polar decomposition = bJ{x)\TYjr{^x)\j where U{x) : C” — Jy\){x) C is 

(the unique) partial isometry with ker[/(x) = kerTpjp^^:) for a.e. x G T^. Hence, in this case 
U{x) = Trj'(a;) {G^/‘^{x)Y and therefore [/(•) : ^ L(C"',f’^(Z^)) is a well defined measurable field 

of partial isometries. Then, Vj : f’^(Z^) given by Vj{x) = U{x) Uj{x) G JyviY for j G In and 

X G are measurable vector fields such that {uj(x)}jg][^ is an orthonormal system in Jw(x), for 
a.e. X G T^; moreover, [5'p;(_7r)]3; Vj{x) = Xj(x) Vj(x) for j G Ir and a.e. x G T^. Since rk[5p;(_7r)]a; = r 
for a.e. x G T^, then we see that Eq. (USD holds in this case. 

Finally, notice that if E{E) is a frame for W then we should have that r = rk[5p;(j-)]j, = dim Jyy(x) 
for a.e. x G which shows the last part of the statement. □ 

Remark 3.3. Let E = be a finite sequence in L^(R^) such that E[F) is a Bessel sequence. 

By Lemma 13.21 there exist measurable vectors fields Vj : —>■ and measurable functions 
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Xj : —)• R+ such that they verify Eq. (IlSp . In what follows we consider the fine spectral structure 

of E{T) that is the weakly measurable function 

9 X !-)• G (non-increasing sequences) , (19) 

where Aj([S'E(j-)]x) = Xj{x) for j € 1 ^( 2 :) and Aj([S'E(jr)]a; = 0 for j > r{x) -|- 1, for x G T^. Hence, 
(Aj([<S'E(j-)]x) )jGN coincides with the sequence of eigenvalues of the positive semidefinite finite rank 
operator = <S'r^(a;) G L(£^(Z^)), counting multiplicities and arranged in non-increasing 

order, for a.e. x G T^. A 

Remark 3.4. Consider the notations from Definition 13.11 We now show that the right hand side 
in Eq. (HZD is well defined. Indeed, by Lemma [32] we get a spectral representation of as 

in Eq. (IlSp in terms of the bounded and measurable functions Aj(-) : —>■ M"*", for j G In- If we 

consider the bounded and measurable function d(x) = dim Jw(x) > r(x} for x G T*' then, using 
Eq. (fT^ we see that 

F^'^^^\rE(x)) = ^ (p{Xj{x)) -|- (d(x) — r(x)) (^(0) for a.e. x G . 


Hence, the non-negative function 


is bounded and measurable and therefore integrable on T^. This shows that the convex potential 
P^{E{IF)) is a well defined non-negative real number. A 


Incidentally, Remark 13.41 above shows that if <^(0) = 0 then the convex potential P^ = P^p does 
not depend on the ESI subspace W. 

Example 3.5. Let W be a ESI subspace of L^(M^) and let P = {/i}iGi„ € If we set (p{x) = x^ 
for X G M"'" then, the corresponding potential on E{P), that we shall denote EP {E{P)), is given by 

FP iE{P)) = f tr(52_^(,)) dx= [ l(r/*(x),r/,(x))|2 dx . 

i,j&n 

Hence, FP {E{P)) is a natural extension of the Benedetto-Fickus frame potential of Eq. (11611 . A 

Remark 3.6. Let W be a SI subspace of L^(]R^) and let A G be a positive operator: in 

[2T] . E. Dutkay introduces the local trace function of A relative to W, denoted rw,yi : [0, 00 ] 

as follows; for x G T^, 

T)/v,a{x) = tr(H [Pyvjx) , 

where tr(-) denotes the usual (semi-finite) trace in L(f'^(Z^)). We can extend the notion of local 
trace function as described above to the following setting: given T G a positive and 

SP operator, we let the local trace function of T with respect to the SI subspace W be given by 

tw,t{x) = tr([r] 3 ; [PW\x) , X gT^ . (20) 

Notice that if H G and T G L{L ‘^is the unique positive and SP operator such 

that \T]x = A for x G then 

xy\;,A{x) = Tyv^rix) , x G . 

If we assume further that W is a ESI subspace, we consider tp G Conv(M''') and take P = G 

W” then 

P^{E{P))=[ rH;,^(5,(^,)(x)dx, 
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where (p{Se{t)) ^ -L(L^(M*^))+ is obtained by the functional calculus. Indeed, notice that in this 
case ip{SE[E)) is a SP operator such that 

MSe{t))]x = [*S'i5(j-)]x) = HSrT{x)) , for a.e. x G . A 

Let W be a FSI subspace. In what follows we show that, under some natural restrictions, the convex 
potentials {E{T)) for finite sequences P G W” detect tight frames for W as their minimizers 
(see Theorem 1,1.91 below!. In turn, this last fact motivates the study of the structure of minimizers 
of convex potentials for finitely generated sequences in (under some restrictions) since these 

minimizers can be considered as natural substitutes of tight frames. In order to state the results 
on this matter, we introduce the following notions and notations. 

Remark 3.7. Let {X, X , fix), (Y, y , fiy) be two measure spaces; we consider their direct sum, 
denoted X 0 T, given by the three-tuple (X 0 T, T 0 T, ® fJ-v), where 

d 

1. X (yY = X UY (the disjoint union of the sets); we further consider the canonical inclusions 
rfx '■ X ^ X ®Y and rfy '■ Y X (BY of X and Y into their disjoint union; hence rjx and 
Tfy are injective functions such that r]x{X) n rjy{Y) = 0 and rix{X) U r]y{Y) = X (BY. 

2. Y = {A® B = r]x{A) U r]y{B) : A € X, B € y}- 

3. fix © fJ-Y is the measure given by fix © fJ^riA ® B) = fix{A) + fiy{B); 

Notice that using the maps rfx and rfy we can consider (as we sometimes do) X,YciX®Y. A 
Notations 3.8. In what follows we consider: 

1. A FSI subspace of L^(M^) of length £, denoted W; 

2. P = G W” such that E{P) is a Bessel sequence; 

3. d{x) = dim Jyv;(x) < £, for x G T^; 

4. The Lebesgue measure on denoted | • | ; Xj = d~^{i) C and pi = |Xj|, i G . 

5. We denote by Cyv = i ■ Pi ■ 

6. The spectrum of W is the measurable set Spec(W’) = Xi = {x G : d(x) / 0}. A 

Theorem 3.9 (Structure of minimizers with norm restrictions). Consider the Notations \S.£^ 
and assume that ll/i|P = 1- If P ^ Conv(M+), then 

P^{E{P))>CyvpiC^^). ( 21 ) 

Moreover, if p ConVs(M+) then equality holds in (f^ iff E{P) is a tight frame for W i.e. 

SEiT) = • ( 22 ) 

Proof. Let {Xij,Yij, \ ■ \ij) where Xij = Xj, Xij = A) the fi-algebra of Lebesgue measurable sets 
in Xj and \ ■ \ij = | ■ |j the Lebesgue measure in Xj, for j G Ij and i B l£. With the notations in 
Remark 13.71 we consider the measure space 

{x,x,fi) = ^^{Xij,Xij, I • \ij) . 

i&e jeh 


12 


For i G and j G I* we further consider the canonical inclusions rjij : Xij —>■ X. Hence, for every 
X £ X there exists unique i G j G Ij and x G Xij = Xi such that r]ij{x) = x. Notice that by 
construction, fj,{X) = i ■ Pi = . 

Let \e{t) • ^ be the measurable function of eigenvalues of E{T) defined as follows: for 

X £ X, let {i, j) G X Ij and x £ Xij = Xi be (uniquely determined) such that r]ij{x) = x; in 
this case we set 

^E{j^){x) = Aj( [5 'e(jp)]£) = Aj(5rjP(i)) , 

where 9 x (Aj([5^(_7r)]a;) )jgis} G is the fine spectral structure of E{X) defined in 

Remark 13.31 We claim that if G Conv(M+), then 

P^{E{E))= [ ^{XEiT){x)) d^i{x). (23) 

Jx 

Indeed, for Eq. (flTl) 

P^{E{E)) = [ P^w(*)(rj'(x)) dx= [ P;!^^^\TP{x)) dx , 

Jt'' Jspec(W) 

where P^^^^\tP{x)) is the convex potential associated with {ip,Jy\){x)) of the finite sequence 
TP{x) = {r fi{x)}i^i^ in t'^(Z^) as defined in Eq. (fT^ (notice that P^'^^^\tP{x)) = 0 for x G 
\ Spec(>V)). Therefore, if x G W for some f G then 

i 

j=i 


For f G If we have that 


f ^ f '^^(X-(SrTix)))dx= [ ip{XE(x)ix))dpix) . 


Therefore, since Spec(VV) = (J X* and X = 0igi^ ©jei; ^i,j ) 

i&c 


P^iEiP)) = Y^ [ P;l^^^\TPix))dx = Y^ [ ^iXE^E)i^))df,ix)= [ ^iXE^E)ix))dpix), 

which proves Eq. (f23]l . In particular, if we take (p{x) = x in Eq. (f23l) we get that 

[ >^E{x)ix) dfj,{x) = f tT{SrE{x))dx= [ '^\\rfi{x)f dx = '^\\fif = 1 . 

Jx Jt>^ Jt^ 


ieln 


Consider the probability measure jl = Cy^ /r. Then, as in Example 12.51 we have that 

/ A£;(;r)(x) (i/2(x) = ■ lx P Xe{E) ( ^ {X, 

Jx 

If we let £ Conv(M'’') then, using the previous facts and Theorem 12.41 we get that 

[2^ 

^ ) = !xJ^ i^w ■ ^x) dp < (p(Xe(e)(x}) dp(x) 


( 24 ) 


= fx <p(Me)(^)) dp(x) # P^{E{P)) , 
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which proves Eq. (|2T|) . If y? E ConVs(K+) and also {E{F)) = ip{Cy^) Cy\; , using Eq. (|23|) . we 
get that 

/ dfi{x) = / ip{C:^^) dfi. 

Jx Jx 

Hence, by Proposition 12.61 and the majorization relation in Eq. (1241) . 

i^EiX))* = Cw 1[0,1] ^ for i E and a.e. x E . 

Therefore, Se(e) = i-e. E{P) is a tight frame for W. Conversely, notice that if Se(e) = 

Py\) then lower bound in Eq. (12X1) is attained. □ 

4 Fine spectral structure of shift generated oblique duals in FSI 
subspaces 

Throughout this section V, W C denote FSI subspaces such that V 0 W-*“ = and 

F = {fi}i£in £ denotes a finite sequence such that E{F) is a frame for W. 

Next we recall some characterizations of the condition S 0 T"*" = for SI subspaces and a 

characterization of shift generated (SG) oblique duals of E{F); these results together with [5] allow 
us to obtain the exact value of the aliasing norm corresponding to the consistent sampling induced 
by the FSI subspaces V and W. In Section 14.21 we obtain a detailed description of the fine spectral 
structure (i.e. eigenvalues) of the frame operators of SG oblique V-duals of the (fixed) frame E{F) 
for W. We will apply these results in Section [SJ where we compute SG oblique dual frames with 
norm restrictions that simultaneously minimize the convex potentials P^ for all (p E Gonv(M'’'). 

4.1 SG oblique duals and aliasing in FSI snbspaces 

Following [25] (see also |T7l [251ES] ) we consider the set of SG V-duals of E{F): 

V^^{F) = V^^{E{F)) = {E{g) E Vv{E{F)) : gl = E V’^} . (25) 

In case V = W then we write 'D^^{F) = {E{F)) (which is the class of SG duals of type I, in the 

terminology of [28]). Given E{Q) E Py^{F) we obtain the following (structured) reconstruction 
formulas: for every / E W and 5 E V, 

/= {f,Tegi)Tifi and 9 = Y. id, Te fi) Te g^. 

(£,i)GZ'=xI„ 

Next we describe some results related with the general assumption for studying oblique duality, 
namely V0 W"*" = L^(M^), for the FSI subspaces V and W, as well as SG oblique duality. The next 
two results can be derived using combinations of results and techniques in [21 [521(55] . 

Lemma 4.1. With the previous notations and assumptions, let Jy and Jyy denote the range func¬ 
tions of the SI subspaces V and W, respectively. Then, 

1. Vy-*- is a SI subspace with range function Jy\;±{x) = [Jvy(x)]-*- for a.e. x E T^; 

2. If Q = Py ^/yyx then Q is a shift preserving operator; 

3. Jp{x) © Jw{x)^ = and [Q]^ = Pj^{x)//Jw{xp for a.e. x E 

4. E(g) E V^^(F) ^ rg(x) is Mx) - dual of TF{x), for a.e xeT^. □ 
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Remark 4.2. Let S and T be closed subspaces of In order to characterize when the 

(algebraic) sum of these subspaces is a closed subspace we recall the Dixmier angle between S and 
T, denoted by [5, T]d € [OjTt], given by 


cos[5,7^d = sup{|(u, w)\,v e Si,w eTi} , 


( 26 ) 


where 5i = {/ G 5 : ||/|| = 1} (and similar for 71). It is well known (see [20]) that [S,T]d > 0 if 
and only if <S Pi T = {0} and <S + T is a closed subspace of 


Assume further that <S © T 
Then (see [20] ) 


L^(R^) and let <3 = TIs//r be tbe corresponding oblique projection. 


IIQII 


1 

sin[5,7lD 


A 


Proposition 4.3. Let S,T © L^(]R^) he SI subspaces of The following statements are 

equivalent: 


1 . 5©r^ = L2(R^); 

2. Js{x) © Jt{x)^ = £^(Z^) for a.e. x G and esssupa-g^fe||.Pj 5 (x)//Jr(x)-LII < oo/ 

3. Js{x)-^ n Jt{x) = {0} and essinf 2 ,gTpfc [Js{x),Jr{x)-^]D > 0. 

In this case we have that [S,T'^]d = ess inf 3 ,gTfe [7<s(a;), Jt{x)'^]d ■ D 


As an application of the previous results we compute the exact value of the aliasing norm (see 
[22[ l30] ) in terms of the relative geometry of the FSI subspaces V and W. Indeed, recall that the 
aliasing norm corresponding to the consistent sampling 

f ^ f = Py^i(vxf : for / G L^(R^) 

denoted A(V, W), is given by 

A(V, W) = sup P^xll. (27) 

eeW-L ll^ll 

The aliasing norm is a measure of the incidence of W'^ in the consistent sampling induced by 
Pyy/and it plays a role in applications of oblique duality. 

Definition 4.4. Let S, T C L^(R^) be closed subspaces. We define the aperture between S and 
T, denoted [5,7^“ G [0,7r/2], as the angle given by 

Remark 4.5. With the notations of Definition [TT] we point out that the aperture [5, T]°' coincides 
with the notion of angle between the subspaces <S and T as dehned in [JH] (and cos([<S,T]“) is also 
known as the infimum cosine angle from [32] ). It is known that the following relation holds (see 

[SaEH]): 

cos([5,7li?)2 = l-cos([5,r]“)2 ^ [5,7l“ = 7r/2-[5,r^]fl. 

Hence, using the relations above and Proposition 14.31 (see also [32]) we get that if <S, T C L^(R*^) 
are SI subspaces such that S © T'^ = L^(R^) then 

[5,7^“ = esssup^fzjk[Js{x),Jrix)]'^ < tiI2. A 
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Consider again the notations of Definition 14.41 and assume further that = 5 © T"*". Then, 

using Remarks 14.21 and 14.51 we see that 


gin[5,r^]B cos[5,7l“‘ 

From this we obtain the following upper bound for the aliasing norm (see m) 

A{S, T) < \\Ps//tA\ = • 

Notice that this known bound is not sharp; indeed, if we take S = T then ^(5, T) = 0 but 
cos[<S,7l“ = 1. 

Next we compute the exact value of the aliasing norm. 

Proposition 4.6. With the previous notations and assumptions, the aliasing norm A(y, W) cor¬ 
responding to the FSI oblique pair (V, W) is given by 

A{V, W) = tan([V, W]“). 

Proof. Notice that by assumption J\;{x) and J>v(x) are finite dimensional subspaces of ) and, 
by Proposition 14.31 we see that J^{x)^ © Jy\){x) = £^(Z^), for a.e. x G T^. Hence, we can apply 
the results from [5], and conclude that 

A{Jv{x), Jw{x)) = \\Pjy,{x)//Mxp Pjyvixp II = tan([Jv(x), Jyv{x)f) , for a.e. x G T''. 
Therefore, using Remark 14.51 we get that 

^(V, W) = ||Pvv;//V-L-fvv-L II = esssup^-g-j-fc tan([Jv(x), Jvv;(a;)]“) = tan([V, W]“). □ 

Conjecture 4.7. We conjecture that Proposition 14.61 holds for the consistent sampling correspond¬ 
ing to an oblique decomposition S © T"*" = in an arbitrary Hilbert space PL. By the results from 
[5] the conjecture holds for finite dimensional S and T. By Proposition 14.61 this conjecture holds 
for some infinite dimensional subspaces S and T as well. A 

4.2 Fine spectral structure of SG oblique duals 

Let E{Q) G and let S'_E(g) denote the frame operator of E{Q). Recall that in this case SE(g) 

is a shift preserving (SP) operator such that [SE{g)]x = ‘S'rg(a;) for a.e. x G and the fine spectral 

structure of E{Q) is the function 3 x ("^i([‘S'_E( 5 )]a:) )jGN € that consists of the 

sequence of eigenvalues of the positive finite rank operator [SE(g)]x = Srg(x)i counting multiplicities 
and arranged in non-increasing order for a.e. x G (see Remark 13.31) . 

In the next result we consider the measurable function d : —>■ {0,... , n} given by d{x) = 

dim Jw(x) = dim Jv(x) for a.e. x G T^. 

Lemma 4.8. Let Q = G be such that E{Q) is a frame for V. Let B G L(L^(]R*^))+ be 

a shift preserving operator sueh that R{B) C V. Then, there exists Z = G sueh that 

B = Se(z) o,nd TE(g) Pe{Z) ~ ^ '^f */rk([B]a;) < n — d{x) for a.e. x G T^. 

Proof. First notice that by considering a convenient finite partition of into measurable sets we 
can assume, without loss of generality, that d{x) = d G N for a.e. x G . Notice that in this case 
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n > d. By Lemma 13.21 there exist measurable vector fields vj : for j G In such that, if 

9 X !-)• (Aj(x) £ (-^+(2^) )'*' denotes the fine spectral structure of E{Q), then 

i^E{g)]x = Y1 Xj{x) Vj{x) ® Vj{x) , for a.e. x G . (28) 

j&d 

Moreover, in this case {vj{x)}j^i^ is an ONB of Jy{x) for a.e. x G T^. Assume that B G 
is a shift preserving operator such that R{B) C V and such that rk([B]a;) <n — d for a.e. x G T^. 
Since [B]x G L(£^(Z^))+ is such that R{[B]x) C Jv(x) for a.e. x G then, using the measurable 
vector fields {vj}j^i^ as above (indeed, the measurable field of matrix representations of [B]^ with 
respect to {vj{x)}j^i^ and the results from [l3]) we get measurable fields Wj : —>■ £^(Z^) for 

j G Id, such that {'Wj{x)}j,=i^ is an ONB of Jv{x) and [B\xWj{x) = Xj{{B]x) Wj{x) for j G 1^ and 
a.e. X G In particular, we see that 

rmn{d^n—d} 

= E Xj{[B]x)^^'^ Wj{x) 0 Wj{x) , for a.e. x G . 
i=i 

Consider the measurable field of positive semidefinite matrices Gg : > Aln(C) given by Gg{x) = 

{{Tgi{x), rgj{x)))i^ iein) X G T^. Again by [33], there exist measurable field of vectors Uj : ^ C"" 

for j G In such that for a.e. x G we have that {uj{x)}j(zi^ is an ONB of C”", Gg{x)uj{x) = 

Xj{x) Uj{x) for j G Irf and Gg{x)uj{x) = 0 for d + 1 < j < n (since Gg{x) and [SE{g)]x have the 

same strictly positive eigenvalues). 

Let V :T^ ^ L(C"', f'^(Z^)) be the measurable field of partial isometries given by 

d+1 < j < d + min{d,n-d}, 

( 0 otherwise. 

Hence V{x)V*{x) is the orthogonal projection onto span{r(;j(x) : j G Imin{d,n-(i}} and thus 
[B]x V{x) V*{x) = [B]x for a.e. x G T^; on the other hand, R{V*{x)) = ker B(x)-*- C span{Mj(x) : 
d + 1 < j < n} Trg( 3 ;)H*(x) = 0 for a.e. x G T^. 

For f G In we set Zj G V determined uniquely by r 2 ;j(x) = [B^J‘^y{x)ei for a.e. x G T^, where 

{^i}i&n denotes the canonical ONB of C”. If we set Z = {zi}i^i^ then Tyz{x) = H(x) for 

a.e. X G hence, using Eq. m, we see that 

[TEiG) T*e^z)]x = Trgix)T^z(x) = Trgix)V*ix) [B^% = 0 , for a.e. x G 
On the other hand, notice that 

[SEiz)]x = Srzix) = V{x) V*{x) [B]^^ = [B]^ , for a.e. x G 

Conversely, assume that there exists Z = {zi]i^i^ G V"' such that B = Se(z) and T^gg-^ '^E{z) ~ 
Then, by Eq. (fTO]) . we get that 0 = Ti-g(x) 'Brz{x) hence ^ ^ “ i’k(Trg(a.)) = n — d(x) 

for a.e. x G T^. Therefore, 

rk [Se(z)]x = rk(5r2(a;)) = rk{Trz(x)) = MBrzix)) ^ “ d{x) for a.e. x G . □ 

Definition 4.9. Let Q = be such that E{Q) is a frame for V with frame operator A = SE{g)- 

Recall that d(x) = dim Jv(x) for x G Then, we consider 

Uv{E{g )) = |a + B : R G L{L^{R'^))+ is SP, R{B) C V, rk([B]a,) <n- d(x), for a.e. x G }. 

A 
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Proposition 4.10. Let E{T)y = denote the canonical V-dual of T. Then, 

{SEis) : E{g) € V^^{E)} = Uv{E{F)*) . (29) 

Proof. Let G = {gi}i£i„ € V” be such that E{g) G V^^{E). Let Z = G V” be given 

II 

by Zi = Qi — /y ^ for i G In ■ Then L'(.H) = {7}i)ez''xn„ is a Bessel sequence in V such that 
r£;(g) = T^^j,.^#+Te(z)- Ill this case we have that Tb( 2 ) T^(jr) = 0 and therefore Tb(^) '^e{t)* ^ 

since Thus, 

SE(g) = i^EiT)* + '^E(Z)) i^EiT)* + = ^E{T)* + ‘^£^(• 2 ) ' 

We conclude that B = Se{z) S L(L^(M^))+ is SP, R{Se(z)) C V and, by Lemma 14.81 that 
rk [Se(z)]x < n — d{x) for a.e. x G T^. 

Conversely, if S' G U^{E{Ff^) then S = + B, where B G L{L‘^(M.^))~^ is SP, R{B) C V and 

rk([B]a;) < n — d{x) for a.e. x € X. By Lemma ITS! we see that there exists Z = such that 

Te{z) 'd^EiT) = ii B = Se{z)- If we let G = be given by 5 ^ = /^ • + Zi for i G then 

E{G) is a Bessel sequence in V such that Tei^q) = ETeiz)- Using that Te(z) = 0 we 

conclude, as before, that 

^E{g) = ^E{J^)* “b ^E{Z) = ^E{J^)* E R = ^ lU 

Proposition 14.101 shows that the set of frame operators of SG V-duals of a fixed frame E can be 
described in terms of the additive model Uy{E{E)^) introduced in Definition 14.91 It turns out 

II 

that the fine spectral structure of the elements of U\;{E{E)^) can be described using a natural 
extension of the Fan-Pall interlacing theorem for measurable fields of positive matrices. We develop 
both results in the Appendix section (see Theorems 16.31 and I6.4p . As a consequence we obtain the 
following 

Theorem 4.11 (Fine spectral structure of V-duals). Let E{E)y be the canonical V-dual frame of 
E{E). Denote the fine spectral structure of E{E)^ by 3 x {X^ ■{x))i^^, x G T^. Let m 
be the measurable function given by m{x) = 2d{x) — n, for x G T^. Given a measurable function 
p : —>■ (£^(N)+)I (decreasing sequences) described as g, = {gi)i(z^, the following are equivalent: 

1. There exists E{G) G 'D^^{E) such that g{x) = A( [SE(g)]x) = A(S'rg(j,)) for every x G T^. 

2. For a.e. x ^ Spec{V), g{x) = 0. For a.e. x G SpeciV), gi{x) = 0 for i > d{x) + 1 and 

(a) in case that m{x) < 0, then gi{x) X Ay j(x) for i G ^d(x)'! 

II 

(b) in case that m{x) > 1, then gi{x) X Ay j(x) for i G 1 ^( 3 ,) and 

hn—d{x)+iiT) Td{x)—m{x)+i(.^') — ^V, ^ ^ ^m{x) • 

Proof. It follows from Proposition 14.IOI and Theorem 16.41 □ 

As a consequence of the description of the fine spectral structure of V-duals of E{F) we characterize 
the existence of tight V-duals of E{F) that are shift generated (compare with [25]). 

Corollary 4.12. With the notations of Theorem \4.11\ then there exists a c-tight V-dual E{G) G 
Vy^{E) if and only if 

d- ^E{T)* ^ ^ 
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2. rk( [c- Pv — < min{(i(x), n — d(x)} for a.e. x € Spec{V). 

Proof. Theorem 14.111 imply that there exists a V-dual E{Q) € Vy^{P) such that SE{g) = c - Py \i 

ff ij 

and only if c > Ay ^(x) for i € and Ay ^{x) = c for i € Ilm(x) whenever m(x) = 2 d{x) — n > 1, 
for a.e. x G Spec(V). These last two conditions are equivalent to the fact that c- P\; > and 

)y 

rk( [c- P\; — < d{x) — m{x) = n — d{x) whenever m(x) > 1 . 

Also notice that in case m(x) < 0 then n — d{x) > d{x) = dim Jv(x). The proof follows from these 
remarks. □ 


Remark 4.13. Consider the notations of Theorem 14.111 As a consequence of Corollary 14.121 we 
get the following dichotomy related with the existence of tight oblique V-duals of E{P): 

1. If n > 2d{x) for a.e. x G T*' then for every c > there exists E{Q) G V^^{P) that is 

a c-tight frame for V. 

2. If there exists C with positive Lebesgue measure such that n < 2d{x) for a.e. x G N 

and there exists a c-tight frame E{G) G Py^(T') then c = A 


5 Applications: optimal oblique SG-duals with norm restrictions 

As before, we consider two FSI subspaces V and W such that >V-*-©V = and P = {/i}iGn„ G 

W” such that E{P) is a frame for W. 

As a consequence of the description of the fine spectral structure of elements in we see that 

the canonical V-dual is optimal with respect to several criteria. Nevertheless, in applied situations, 
the canonical dual might not be the best choice: for example, we can be interested in duals of E[P) 
such that the spectrum of their frame operators are as concentrated as possible. Ideally, we would 
search for tight dual frames for E{P), although Corollary 14.121 shows that there are restrictions for 
the existence of such duals. 

In order to search for alternate V-duals that are spectrally more stable, we proceed as follows: for 
w > EiGin where E{P)* = we consider 

V^%{P)=V^%iE{P)) {E{g)GV^^{P): G = and > W ■ 

iein 


Notice that if rc > EiGin then E{P)y ^ and therefore, it is natural to ask whether 

there is an optimal dual fulfilling the previous requirements. Using the identity 


i&n 


Y1 dx 

i&n 


[ tr( [SE(g)]x) 
XT* 



'^fii{x) dx 

iGN 


(30) 


where A( [^^(g)]®) = for a.e. x G we see that Theorem 14.111 gives a complete 

solution to a frame design problem in the sense that it allows to get a complete description of the 
eigenvalue lists of the frame operators of elements in It is then natural to seek for those 

oblique SG-duals E{Q) G Py‘E(-^) that minimize the convex potentials P^, for ip G Conv(M+); 
in order to deal with this problem we first examine a construction known as water-filling in terms 
of submajorization, in the general context of measure spaces (see Theorem 15.51) . We then apply 
these results together with the properties of submajorization and results from matrix analysis to 
conclude that there are structural optimal duals with norm restrictions. These optimal solutions 
are obtained in terms of a non-commutative water-filling construction. 
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5.1 Water-filling in measure spaces 

The water-filling construction goes back to the work of Shanon [l6], as the solution of an optimal 
spectral allocation problem (see m)- The water-filling strategy has also been the main tool in the 
design of channels with optimal capacity (see m and the more recent work on iterative water-filling 

techniques mm)- 

As a first step towards an extension of this construction, we examine its scalar counter-part in 
the general context of measure spaces. In the next section we show that the water-filling technique 
produces optimal solutions in the general (non-commutative) context of measurable fields of positive 
semidefinite matrix valued functions. 

Throughout this section the triple {X, X, fi) denotes a probability space. Recall that we denote by 
L-(A,/r)+ = {/GL-(A,^) :/>0}. 

Definition 5.1 (Water-filling at level c). Let / G L°°(A,Given c > essinf/ > 0 we consider 
fc G T°°(A, /r)"'' given by fc = max{/, c} = / -|- (c — f)^, where denotes the positive part of a 
real function g. A 


In order to study the submajorization properties of the function fc obtained by the water-filling 
construction as above, we consider the following result in which we obtain a simple relation between 
the decreasing rearrangements of / and fc- 

Lemma 5.2. Let f G L°°{X,g,)^ and let c > essinf/ > 0. Consider the number 


So = /.i{x G X : f{x) > c} . Then fc{s) 


f*{s) if 0<s<so; 
c if So < 'S < 1 • 


Proof. Notice that by Eq. m, for 0 < s < So we have that 


ns) 


sup {t G M+ 
sup {t G M+ 
sup {t G M+ 
sup {t G M+ 


g-{x G X : f{x) > t} > s} 

g,{x G X : f{x) > t} > s and t > c} 

/i{x G X : fc{x) > t} > s and t > c} 

fi{x G A : fc{x) >t}> s} = fc{s) . 


(31) 


It is straightforward to see that if so < s < 1 then //(s) = c. 


□ 


In order to prove Theorem 15.51 below, we shall need an explicit statement of some re-parametrized 
versions of the basics results of section 12.31 

Lemma 5.3. Let a, b gR be such that a < b and let k G L°°{[a ,6], i')~^ be a non-increasing right 
continuous function, where u = (b — a)~^ dt is the normalized Lebesgue measure on [a, b\. Then 

1. The decreasing rearrangement k*(t) = k(^ {b — a) t a) for every t G [0 , 1). 

2. Fix a constant c G M. Then 



for every s G [a, 6] . 


Proof. Straightforward. □ 

With the notations of Lemma 15.31 above, notice that item 2. is a restatement (using the re- 
parametrization from item 1) of the submajorization inequalities corresponding to c -<w k in 
([a, 6], n) whenever c < k du (see Example 12.51) . 
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Remark 5.4. Let / E and consider (p : [ess inf/, oo) ^ given by 


He) = [ fcdn= [ f{x) + (c - /(x))+ dnix) 
Jx J X 


Then, it is easy to see that cp has the following properties: 

1. (/(ess inf/) = and linic->.+oo </(c) = Too; 

2. (p is continuous and strictly increasing. 

Hence, for every w > f dfi there exists a unique c{w) = c> ess inf/ such that 


(p{c{w)) = w i.e. / fc(w) dfi = w 
Jx 


(32) 


A 


Theorem 5.5 (^^-optimality of water-filling). Let f E , take w > f dfi and consider 

the constant c{w) = c as in Remark\5.4\ Then, for every h E L°°(X,//)+, 


f < h and 


h d^> w 


fc h . 


IX 


Proof. Assume that f < h and h dfj.> w. If we let sq = ^ ^ • fH) > H then, by Lemma 

15.21 we have that Eq. (1311) holds. Thus, using Remark 12.21 we see that if 0 < s < sq then 


ffft) dt= f*it) dt< h*{t) dt . 


(33) 


Fix now So < s < 1 and consider 

/■So rso 1 

a = h*{t) dt — /c (^) dt > 0 and k = h* + - -a E L°°([0,1], dt)~ 

Jo Jo 1 - So 

Notice that A: is a non-increasing right continuous map. In this case we get that 

f k{t) dt = f h*{t) dt + a 
J So J Sq 

/■I / rso rso ^ 

= / h*{t)dt+[ h*{t)dt- feit) dt 

J so \Jo Jo J 


> tc - (rc - (1 - so)c) = (1 - so)c . 

Then, Lemma E2] (applied to the map A:|[gg q ) implies that 

(s — So) c < f k{t) dt = f h*{t) dt + - -^ a for every s E [sq , 1] 

Jso Jso 1 - So 

Hence, using the inequality above and Lemma [521 we conclude that for sq < s < 1 

cso 


rs PSQ 

/ ff{t)dt = / h*{t) dt — a + {s — so)c 

Jo Jo 

[ h*{t) dt +(-—^ a < [ h*{t) 

Jo \1 “ So / Jo 


< 

/o ' " VI “ So 

This last fact together with Eq. ([55|1 show that fc -<w h. 


dt . 


□ 
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Theorem [53] above implies a family of integral inequalities in terms of convex functions involving 
the water-filling of a function / at level c. We will need these facts in order to show the optimality 
properties of the non-commutative version of waterfilling. 

Corollary 5.6. With the notations of Theorem \5.51 if (f G Conv(M+) is non-decreasing then 


/ ip o h dgL> / if o fc dp . 

Jx Jx 


(34) 


If there is a non-decreasing ip G ConVs(M"'') such that equality holds in Eq. (l33l then h = fc- 


Proof. The first claim is a consequence of the submajorization relation in Theorem [53] and Theorem 
12.41 If we further assume that ip € ConVs(M+) is such that equality holds in Eq. (f34|l then, by 
Proposition 12.61 we see that /* = h*. Let B = {x G X : f{x) > c} so that sq = p{B). Then, it is 
straightforward to show that 


{f-lBYis) 


f*{s) if s€[0,so); 
0 if s € [so, 1). 


Notice that, in particular, (/-Is)* = l[o,so)'/c' On the other hand we have {hTs)* = l[o, 
Hence, since h>h-lB>fXB,hy Remark 12.21 we have that 


h*>{h- 1b)* = l[o,.o) ■ {h ■ 1b)* > (/ • 1b)* = l[o,so) ■ fc (h ■ 1b)* = (/ • 1b)* . 

Therefore, again by Remark 12.21 we get that h ■ 1 b = f ■ 1 b = fc Xb > c ■ 1 b, where the last facts 
follow from Definition 15.11 Finally, notice that 

T{h-\{c})) = \{h*)-\{c})\ = \{f:)-\{c})\ = 1 - p{B) , 

which shows that h ■ 1x\b = c • lx\B and hence h = fc- □ 


5.2 Optimal SG-duals with norm restrictions: NC water-filling 

In what follows we show the existence of structural optimal SG oblique duals of a fixed frame E{IF) 
with norm restrictions, as described at the beginning of Section [5] That is, we explicitly construct 
a dual frame E{Q°'^) G such that for every E{Q) G T’y‘^^(3) then 

p^(E(g-p))<p^(E(g)), 

for every convex potential associated to a non-decreasing ip G Conv(]R'*’). Moreover, the ar¬ 
guments involved in this construction show that (structural) optimal SG oblique duals with norm 
restrictions share several spectral properties. We end the section with a non-commutative counter¬ 
part of the water-filling construction for functions, that allows to describe the spectral and geomet¬ 
rical structure of optimal SG oblique duals in Dy^^(3) in some detail. 

Theorem 5.7 (Optimal duals in T>y‘^^(3)). Let V and W be ESI subspaces of such that 

W-*- © V = L^(]R^). Let T = be such that E{IF) is a frame for W and w > 

where E{P)y = {Tf Then, there exists G V” such that: 

1. E{g-P) G and 

2. For every g = {gi}iei„ such that E{g) G and every non-decreasing ip G Gonv(M+) 

we have that 

p'^{E{g^n) < p^{E{g)). 
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Proof. Let d{x) = dimJv(x) = dimJw(x) for x € T^. For each i £ In, let Xi = d~^{i) C T^, 
Pi = \Xi\ (the Lebesgue measure of Xi) and rj = min{n — i, i}. Since E{IF) is a frame for W 
then Spec(V) = Spec(W) = Also, for i £ In and j £ I,., we consider the measure space 

{Xij,Xij, I • \ij), where Xij = Xi, Xij = Xi denotes the fi-algebra of Lebesgue measurable sets in Xi 
and I ■ = I • \i denotes the Lebesgue measure in Xi. Then, using Remark I8.7I we construct the 

measure space 

{Y,y,u) = ^^{x,j,yij,\-\,j). 

In particular, I'fV) = r* -pi. We further consider the canonical inclusion maps pij : Xij —>■ Y 

for i £ In and j £ ■ 

Let Q = be such that E{Q) £ We shall denote by A = and S = SE(g) € 

L(L^(M^))+ . By Proposition 14.101 S = S' , + B = A + B, loi some B £ L(L^(M^))+ which 

is SP, R{B) C V and rk([i?]a;) < n — d{x) for a.e. x £ T^. Let i £ I„; using Lidskii’s additive 
inequality (see 0 ) we get that for a.e. x £ Xi 

()jgj. ^ (-^j( ['S']x) )jgj. , (35) 

while Aj( [S]a;) = 0 for j > i + 1. Notice that R{[B]x) C Jv{x) and rk([i?] 3 ;) < n — i. Therefore 
rk([B]a;) < min{n — i, i} = ri. Then, for x £ Xi we have that 


Ai-j + l( [A]^;) + Xj{[B]x) 


Ai_j+i([A]3,) + Xj{[B]x) if 1 <j <ri; 
Ai_j+i( [AJa;) if ri + l<j<i. 


(36) 


Now, Eq. (1351) together with Eq. p6]l imply that, for any ip £ Conv(R+): for a.e. x £ Xi then 


Ti i 

^V^(Ai-j + l( [A]a;) + Aj([B]3;)) + ^ y5(Ai_j + l( [A]a;) ) < ^ </?( Aj([S]a;) ) . 


i=i 


j=ri+l 


j&i 


(37) 


With the previous notations, we now consider the measurable function h : Y ^ defined as 
follows: for x £ Y, let (i, j) € In x 1^. and x £ Xij = Xi be (uniquely determined) such that 
Pij{x) = x; in this case we set h{x) = Ai_j+i( [A]^) + Xj{[B]x). If we let wq = 
we assume that E(Q) £ T’y‘^^(J^) then, using Eq. (f30]l we see that 



tT{[B]x) dx 



— [A]a;) dx > w 


Wq > 0 . 


Consider now the measurable function f : Y ^ M+ given by f{x) = Ai_j+i( [AJ^) for x £ Xij = Xi, 
with (i, j) € In X Ir^ such that Pij{x) = x. Arguing as in the proof of Theorem 13.91 we get that 



di^ = '^ f ^Ai_j+i([A] 

iein ieir, 


;) dx . 


(38) 


Moreover, the previous facts show that if E{Q) £ Vy^^(F) we have that h > f and 



Wo) + 



f{x) dv 


def / 

= w . 


Let c = c{jj^) be as in Remark 15.41 and consider fc as in Definition 15.11 both with respect to the 
probability space (Y, y, v), where v = v{Y)~^ v. By Corollary 15.61 and the previous remarks we see 
that \i ip £ Conv(]R'’') is non-decreasing then 


/a 


ip o dv < / ip o h dv . 


lY 


(39) 
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For j € In we consider the measurable functions : Spec(V) —)■ M"'' defined as follows: for i G In 
and X G Xi, 


^j(^) 


fciVijix)) = max{c, Xi-j+l{ [^]a:)} 

[-^x) 

0 


if 1 < J < r* 
if Xi + 1 < j < i 
if * + 1 < j < n 


Notice that by construction, if (/? G Conv(M'’') then 



dx. 


(40) 


(41) 


Using the definition of / and the properties of fc from Remark 15.41 we see that if x € Xj, then 
there exist A°’’(x) > ... > A°f(x) > 0 such that 


(j{x) = Xi-j+ii [^]x) + A°P(x) for every j < Xi . (42) 

Let —)> .^^(N)+ such that: fij{x) = 0 for j € N whenever x € \ Spec(V), while 

for i € In and x € Xi then fij{x) = 0 for j > i + 1 and 

{d'jix))j&i = [(0(a;))ieij'‘'- (43) 


Putting the previous remarks together we see that fi = (^j)jgH satisfies the conditions of item 2 in 
Theorem 14.Ill Thus, there exists = {5'°’’}iGin such that E(Q°p) G and A( [^^(gop)]^;) = 

(/Uj(x))jgis} for a.e. x G T^. In this case, if we consider Eq. (I30l) . use Eqs. (|i2]l . (|43]) and we take 
(p{x) = X in Eq. (|ini we have that 



where we have also used the relation in Eq. ()38p above. In particular, satisfies item 1. in the 
statement. Now, if E{Q) G using Eqs. 1(371) . (f39]l . (fl0|) and (|43]l then, 


P^{E{Q)) > / ipohdi^ + E E (p o Xi_j_^i(^^A.^x') dx 

3=ri+l 

> f ^°fcdu + '^f ^ Xi-j+i{[A]^) dx = P^{E{g°P)) 

ieln j=ri+l 

where we have also used Eq. m and the fact that A( [(*S'£(gop)]a;) = /u(x) for a.e. x G T^. □ 

Corollary 5.8 (Essential uniqueness of optimal V-duals with norm restrictions). With the notations 
of Theoxem [13 and its pxoof, assume that Q = is such that E{g) G Vy^^{E) and that thexe 

exists a non-decxeasing ip G ConVs(M''') such that 

P^{E{g^'^)) = P^{E{G)). (44) 

Let B G L{LP‘(R^))^ he SP, with R{B) C V and such that Se^q) = ^eF)* + -^ = ^ + R. Then, 
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EiGi„ IlftiP = 

2. There exist c > 0 and measurable vector fields Uj : —>■ for i ^ In such that 

is an ONE of Jv{x) for a.e. x € Spec{V), 

[5'g(jp)#]x = [A]x = ^ Ai([^]j:) Vi{x) i^vfix) , for a.e. x € SpeciV) 
and such that for a.e. x E Spec{V) we have that 

d{x) 

[B]x= ^ {c-Xi{[A]x))'^ Vi{x) ^Vi{x) , 

i=r(x)+l 


where r{x) = max{2(i(x) — n, 0}, for x € T^. 

The constant c (= > 0 does not depend on Q. Moreover, in this case P^{E{Q)) = 

{E{Q°'^)) for every non-decreasing if € Conv(M"*“). 

Proof. We use the notions and notations from the proof of Theorem 15.71 Arguing as in the last 
part of the proof of Theorem 15.71 we see that Eq. (1441) implies that 

/ If o h du = / tp o f^ di>, 

Jy Jy 

where u = h'{Y)~^ u is the probability measure obtained by normalization of n. By Corollary 15.61 
we get that h = fc, where c = is as in the proof of Theorem 15.71 Therefore, for i € In and 

X G Xi, then 


\ '1 , \ (\R] \ - / max{Ai_j+i([A],r), c} if 1 < j ; 

+ XmV - I a._,^i([A].) if r, + 1 < j < z. 

Moreover, by Eq. (EZD and the properties of h we have that 

P^iE{G)) = E/ 

iGln j&i 

> f ip oh dv + oX,.j+ii[A]x)dx = P^{Eig)) . 

„v-Tr d X.i ■ —I 1 


(45) 


• j=ri+l 


Therefore, we should have equality Eq. (I3ZD for a.e. x G Xi and i Gin. Since ip is strictly convex, 
then the majorization relation in Eq. (1351) together with the case of equality in Lidskii’s inequality 
(see the Appendix section in [39]) imply that for i G In and a.e. x G Xi there exists an ONB 
{zj{x)}j£l^ of Jv{x) (but not necessarily of measurable vector fields as functions of x) such that 


= ^Aj([A],r) Zj{x) ®Zj{x) and [B]x = Xi-j+i{[B]x) Zj(x) ® Zj{x) . (46) 

j&i j&i 


Let P G L(L^(M^))+ denote the orthogonal projection onto TZ = R{B), so that P is SP and [P]x = 
Pr{IB]x) every x € Let p : —)> N be the measurable function given by p{x) = tr( [Pja;) for 
X G T^. Then, by inspection of Eqs. (1451) and ([16]) we see that P P, 

[P]x [Al]. [P]x + [B]x = C • [P]x , for a.e. x G 
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and, for i and x € Xj, 

i—p(x) 

[I - P]x [A]^ [I - P]x= ^ ^ji[A]x) Zj{x) (g> Zj{x). 

i=i 

Since o + (c — a)~^ = maxja, c} for a, c > 0, these last facts imply the existence of measurable 

vector fields Uj : —>■ for i G In with the desired properties; indeed, the previous identities 

show that we just have to consider measurable fields of eigenvectors of the operators P 

_ )\? 

and (I — P) S'n,/T-N# (I — P), whose existence follow from Lemma 13.21 

Finally, if is as in Theorem 15.71 a careful inspection of the proof of that theorem shows that 
K [SE{g)]x) = A( [^^(gop)]^;) , for a.e. x G SpecCF), 
which implies the optimality properties of E{G) for a non-decreasing ^ G Conv(M+). □ 

Notice that with the notations of Corollary 15.81 above, we see that for a.e. x G then 

r{x) d{x) 

[SE(g°p)]x = E Aj([^] 3 ;) Vi{x) (g) Vi{x) + ^ max{Aj([^] 3 ;), c} Vi{x) (g) Vi(x ), 

i=l i=r(x)+l 

where we have used that a -|- (c — a)"*" = max{a, c} for a, c > 0. In particular, notice that 

^d{x)[SE(g°p)]x P max{c, ([^j^)} 

which implies that the condition number of [S'£;(gop)]a, is smaller than or equal to the condition 
number of [A]x = - both acting on Jv{x) - for a.e. x G T^. That is, the optimal oblique 

dual E{G°^) improves the (spectral) stability of the canonical oblique dual E{P)y = E{P^). 

The representation of [^^(gop)]^; above motivates the following construction, which also characterizes 
all elements of which are minimal in the sense of Theorem 15.71 

Definition 5.9 (Non-commutative water-filling at level c in U\;{E{G)) )• Let G = {gi}i£i„ be such 
that E{G) is a frame for V with frame operator A = S^ig) ■ By Lemma 13.21 we can consider 
measurable vector fields Vj : for j G In such that 

Mx= E Xji[A]x) Vj{x) (g) Vj{x) (47) 

is a spectral representation of [A]x , where is an ONB of J\;{x) (here d{x) = dim J\;{x) < 

n), for a.e. x G 

Given c > 0 then we define the (non-commutative) water-filling of A at level c with respect to the 
representation in Eq. (113, denoted Ac G U\;{E{G )), as the unique positive SP operator such that 
operator R{Ac) C V and 

d(x) 

[.4c](x) = [Ac]x = ^ \i{[A]x) Vi{x) Vi{x) + ^ max{Ai([.4]j,), c} ni(x) (g)i;i(x) (48) 

ielr(x) i=r{x)+l 

where r(x) = max{2(i(x) — n,,0} (recall Iq = 0) for a.e. x G A 

Remark 5.10. With the notations of Definition 15.91 
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1. We point out that [Ac] as described in Eq. (HSl) is a well defined measurable field of positive 
semidefinite operators that is essentially bounded. 

2. Notice that in the spectral representation of given in Eq. (@8]), the eigenvalues are not 

necessarily arranged in non-increasing order. 

3. Einally notice that Ac G U\;{E{Q)), since [Ac\{x) — [A\x is a positive operator with rank at 

most d{x) — r{x) < n — d{x), for a.e. x G T^. A 

We end this section with the following comments: with the notions and notations of Theorem 15.71 
let A = and consider measurable vector fields tij : —>■ for i G In, such that 

E Ai([^]x) Vi{x) ® Vi{x) (49) 


is a spectral representation of [A]x with respect to an eigen-basis {vi{x)}i^i^^^y where d{x) = 
dim(Jv(a:)), for a.e. x G T^. Let c > 0 be such that, if Ac is the water-filling of A at level c with 
respect to the representation in Eq. (|49p then, 

/ tr([^c](2:)) dx = w . 

By construction Ac G U\;{E{Fy^) and therefore, by Proposition 14.101 there exists Qq G 

such that Sqq = Ac- As we have already noticed, G°'^ from Theorem 15.71 is constructed in this way; 

hence, in this case we have that for every non-decreasing ip G Conv(M“’'), 

P^iEiOo)) < P^{E{g)) , for every gl G . 

Moreover, by Corollary 15.81 any structural optimal frame G G (i-e. such that ^ is a 

P,^-minimizer in for every ip G Conv(M"*“)) is obtained in this way. That is, the structural 

optimal SG V-dual frames for E{E) with norm restrictions are exactly those Q G for which 

their frame operators are obtained in terms of the non-commutative water-filling construction from 
Definition 15.91 

6 Appendix: spectral structure of Uy{E{Q)) 

In what follows we consider a measure space {X, A, p), such that A C is a Lebesgue measurable 
set, X denotes the cr-algebra of Lebesgue sets in X and p is the Lebesgue measure restricted to X. 

Proposition 6.1. Let G ; A —)• be a bounded measurable field of positive semidefinite 

matriees with assoeiated measurable eigenvalues Aj : A —>• M"'' for j G 1^ such that Ai > ... > A„. 
Assume that the measurable funetions : A — >• M"*" for j G I^-i satisfy the interlaeing eonditions 

Xj{x) > fij{x) > Xj+iix) x-a.e. for j G In-i ■ (50) 

Then there exists a measurable map IE : A — >■ A4n,n-i(C) sueh that W*{x) W{x) = In-i and 

A(IE*(x) G(x) IE(x) ) = (/3i(x), (x)) , for a.e. xGX. (51) 

Proof. We argue by induction on n (the size of G). Notice that /?!>...> fin-i by Eq. (150]) . 
Using the results of [33], we can consider measurable vector fields Uj : X ^ for j G 1^ such that 
{uj{x)}j^i^ is an ONB of eigenvectors of G{x) for a.e. x G A. 
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Assume first that fin-i = A^. Set G'{x) = V{x)* G{x) V{x) where V{x) is the n x (n — 1) matrix 
whose columns are the vectors ui{x),... ,Un-i{x), for x G X. Then, G' is a bounded measnrable 
field of (diagonal) positive semidefinite matrices of size n — 1 with measurable eigenvalues Xj : X ^ 
M"*" for j € In-i- If we assume that we can find a measurable function Z : X ^ A4n-i,n-2{C) such 
that Z*{x)Z{x) = In -2 and X{Z*{x)G'(x) Z{x)) = (/3i(x),..., f3n-2ix)) for a.e. x € A, we let 

Wix) = (, for X € A . 

\I'n-2 I / 

Then, it is easy to see that IT : A —)• AI„^„_i(C) has the desired properties. By iterating the 
previous argument and considering a convenient partition of A into measurable sets, we can assume 
without loss of generality that 

Xj{x) > I3j{x) > Aj_|_i(x) , for a.e. x G A , j £ . 


In this case we set 


lj{x) 


n (Ai(x)-/3i(x)) 

iGlfi—i 


n (Aj(x) - Afc(x)) 


for X G A , j £ In ■ 


The previous assumptions (strict interlacing inequalities) imply that 7 j(x) > 0 is defined for a.e. 
X G A; moreover, the functions 7 j : A —>■ M"*" are measurable for j G In- 

Set : X —)• M'*' for j £ I„, let v = ■ X —>■ C"" and let P : A —)■ A4„(C)''“ 

given by P(x) = / — Py(x) (the orthogonal projection onto {n(x)}-'“, notice that v(x) ^ 0 a.e.). Let 
Px{t) G M[t] denote the characteristic polynomial of P{x) G{x) P{x). Then, a well known argnment 
in terms of alternate tensor products (see 0) shows that 

Px{t) = t ^2 n(^ “ Afc(x)) ^ px{Xj{x)) = Xj{x) (Aj(x) - j3i{x)) 

jeln ieln-1 


for a.e. x £ X, j £ In and Px{0) = 0. Therefore, 

Px{t) = t rr ^2 ~ ^ X £ X , 

j&n-l j&n 

by comparing the leading coefficients of the two representations of the polynomial. This last nor¬ 
malization condition shows, in particular, that P(x) = I — v{x) 0 v{x) for x G A a.e. and hence P 
is a measurable function. 

Finally, let {wj : A —>■ C”’}jg][„ be a measurable ONB of eigenvectors functions for P such that 
P{x)wn{x) = 0 for a.e. x G A. Set IT : A —>■ Ad„^„_i(C) such that IT(x) is the n x n — 1 
matrix whose columns are the vectors rci(x),... ,Wn-i{x); then IT is a measurable function with 
the desired properties. □ 

Lemma 6.2. Let Aj : A —>■ M"*" for j £ In he measurable functions such that Xi > ... > Xn. Let 
d £ In_i and let (3j : X ^ M'*' for j G Id be measurable functions such that Pi > . .. > fid and such 
that they satisfy the interlacing inequalities 


Xj{x) > Pj{x) > Xn-d+j{x) , for .a.e. X £ X , j £ld . (52) 

Then, there exist measurable functions 'jij'.X^ M+ for 0 < i < n — d and j £ In-i such that: 

1- 70 ,i = Aj for j £ In and 'jn-dj = Pj for j £ Id; 
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2. For 0 < i < n — d then 'yi,j(x) > 'yij+i(x) for j E In-i-i, for a.e. x E X; 

3. For 0<i<n — d—1 then 'Jijix) > 7 i+i ^j(x) > 7 ^^ j+i(x) for j E In-i-i, for a.e. x E X. 

Proof. We argue by (decreasing) induction in terms of d. Notice that the statement is trivially true 
if d = n — 1. Assume that the result is true for d + 1 interlacing measurable functions for some 
d E In -2 ■ Given the measurable functions j3j for j E 1^ as above, we shall construct measurable 
functions Oj : A —>• M"*" for j £ I(i +1 such that 

Aj > aj > \n-{d+i)+j for j E Id+i and aj > jdj > aj+i for j G Id, (53) 

and hence ai > ... > ad+i ■ Notice that the lemma would be a consequence of this construction 
and the inductive hypothesis (where the maps aj play the role of jn-d+ij for j G Id+i). 

First notice that by the interlacing inequalities in Eq. (|52l) we have that 

min{Ar.+i, /3r} > max{/3r+i , An-d+r} , for r E Id_i and for a.e. x G X . (54) 

for j G Id+i, as follows: 

__ J max{/3j , A„_(rf+i)+j} if 1 < j < d; 


We define aj : X 


a,- := 


min{/3d , Ad+i} if j = d + 1. 


(55) 


By construction the functions aj are measurable, and it is easy to check (by using Eq. ()54l) l that 
they satisfy Eq. (I53l) . □ 

The following result is the Fan-Pall interlacing inequalities theorem for measurable fields of positive 
operators. 

Theorem 6.3. Let G : X ^ A4„(C)'’' he a bounded measurable field of positive semidefinite 
matrices with associated measurable eigenvalues Aj : A —>■ M+ for j G In such that Ai > ... > An- 
Let d G In-i and let fdj : X ^ forj G Id be measurable functions such that /3i > ... > fdd. Then 
the following conditionas are equivalent: 

1. Aj(x) > /3j(x) > \n-d+j{x) , for a.e. x G X , j G Id- 

2. There exists a projection valued measurable function P : X ^ such that 

rkP{x) = d and A(P(x) G(x) P(x)) = (/3i(x),...,/3d(a;), On-d) , for a.e. xGX. 

Proof. Assume first that the functions {/3j}jg]i^ satisfy the interlacing inequalities in item [TJ Let 
'jij : X —)■ M"*" for 0 < i < n — d and j G In-i be measurable functions as in By Lemma 16.21 By 
Proposition 16.11 there exists a measurable function Wi : X —)• Ain, n-i(C) such that 

VFi(x)*VFi(x) = Jn-i and A(kFi(x)*GkFi(x)) = (7i7(x),... , 71 , n-i(a^)) , for a.e. x E A . 

Arguing as before, using Proposition 16.11 we can construct for 2 < i < n — d measurable functions 
Wj : A —> Ain-i+i,n-i{C) such that Wi{x)*Wi{x) = In-i for a-e. x E A and 

X(Wi{x)* • • • kFi(x)*G(x) kFi(x) • • • Wi{x)) = ( 7 i,i(x),... ,'yi^n-iix)) , for a.e. x E A . 

Let W = Wi ■ ■ ■ Wn-d ■ X Aln,d(C) which is measurable by construction and notice that 

W*{x)W{x) = Id and A(kF(x)*G(x) VF(x)) = (/3i(x),...,/3d(x)) , for a.e. x E A . 

Hence, if we set P = WW* : X Ain{C) then P is a measurable field of projections with the 
desired properties. 

Conversely, assume that there exists a projection valued measurable function P : A —>■ Aln(C)^ 
satisfying item [2j Then item [1] is a straightforward consequence of the so-called Cauchy interlacing 
inequalities from matrix analysis (see for example ID). □ 
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Let Q = {(7i}iei„ be such that E{Q) is a frame for the SI subspace V, with frame operator A = SE{g)- 
Recall that (see Definition 14.91) 

Uv{E{g) ) = {A + B:B € L{L‘^{R’^))+ is SP, R{B) C V, rk([R]^) <n- d{x ), for a.e. x € T^} . 

Using the Fan-Pall inequalities for measurable fields of matrices we can now describe the fine 
spectral structure of the elements in U\;{E(g )) 

Theorem 6.4. Let V be a SI subspace in L^(]R^) with Spec{V) C X and let Q = be such 

that E{Q) is a frame for V with frame operator A = SE(g)- Let d : X ^ be the measurable 
function given by d{x) = dim Jv(x), for x € Spec{V), and let m{-) = 2d{-) — n. Given a measurable 
function /x : X —)• (decreasing sequences) the following are equivalent: 

1. There exists C G U^{E{Q )) such that ii{x) = A([^]a;) for a.e. x G X; 

2. g{x) = 0 for every x ^ SpecfV). If x G Spec{V) then fii{x) = 0 for i > d{x) + 1 and 

(a) in case that m{x) < 0, then fii{x) ^ Ai([^] 3 ;) for i G Id{x)i 

(b) in case that m{x) > 1, then gi{x) ^ Ai([^] 3 ;) for i G Id{x) ^.nd 

hn—d{x)+i(.^) Td{x)—m{x)+i(.^) — for i G ^m{x) • 

Proof. First notice that by considering a convenient finite partition of X into measurable sets we 
can assume, without loss of generality, that d{x) = d G N for a.e. x G Spec(V). 

Let C G U\>{E{Q )), and assume that /x = A( [C]). By hypothesis, there exists B G L{L'^(R^))~^ SP, 
with R{B) C V, rk([i?] 3 ;) < n — d{x) for a.e. x G T^, such that C = A + B. By Lemma H78l there 
exists Z = {ziji^i^ G such that TE(g) ^ ^ = Se{z)- If we let ^ -h .2 = {gi + 

then TE{g+z) = '^E{g) + '^E{z) and 

SE{g+z) = TE(g+z) TE{g+z) = ^Eig) + Se{z) = A + Se{z) = C 

with R{[C]x) = Jv{x) and dim Jv(x) = d < n. Then, 

( ['^E{g+z) TE{g+z)]x) = hjix) for j G and a.e. x G Spec(V) . 

Moreover, if we let P : Spec(V) —)■ A4„(C)"'“ be the projection valued measurable function such 
that P{x) is the orthogonal projection onto span {Ft/(x)} = R{T(gf^^.^) then, using again that 
TE[g) '^E{Z) = 0 we see that 

P{x) ( '^E{g+Z)]x ) R(x) = [TE;(g) TE(g)]a; , for a.e. x G Spec(V). 

Since rk(P(x)) = d < n and 

['^E(g) TE(g)]x) = Aj( [TE(g) = >^ji[A]x) , for j G Id and a.e. x G Spec(V) 

then, using Theorem 16.31 we conclude that that the Fan-Pall inequalities hold between 

(Mi(^)) • • • ) Tdix ^) On—d) and (Ai(, ... , Ad(, On—d) ■ 

A careful inspection of these inequalities for the previous vectors shows that the inequalities in 
items 2.a and 2.b. above hold (according to the relation between n and d). 

Conversely, let g, : Spec(V) —>■ £^(N)''' satisfy the conditions in item 2 and let D^(-) : Spec(V) 
A4n(C)+ be the measurable field of positive semidefinite matrices such that D^{x) is the diagonal 
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matrix with main diagonal ..., fJ-dix),0n-d) for x € X. Then, by Theorem 16.31 there exists 

a projection valued measurable function P : Spec(V) —>■ A4niC)~^ such that 

tr(P(x)) = d and X{P{x) D^{x) P{x)) = (Ai([^] 3 ;), ... , Arf([A]j,), On-d) G (M+)"- . 

In this case we see that 

Pix) Dl/^) = (Ai([A],), ... , Ad([^].), 0„_d) € (M+)" . 

Let Dx{x) be the diagonal matrix with main diagonal (Ai([A]a;), ■ ■ ■ , Ad([A] 3 ;), On-d). By taking 
an appropriate measurable field of unitary matrices U{x) : Spec(V) ^ Ad„(C) we conclude that 

Dx{x) = U{x)* P{x) U{x) , for a.e. x G Spec(V). (56) 

Arguing as in the proof of Lemma 14.81 we see that there exist measurable fields of vectors Vj : 
Spec(V) —>■ for j G Id, such that [A\xVj{x) = Aj([A] 3 ;) Vj{x) and B{x) = {vj{x)}i^i^ is an 

ONB of Jv{x) for a.e. x G Spec(V). We finally consider B G S.P. with R{B) C V, 

uniquely determined by the condition: 

= U{xr (Ziy2 (j _ p(^)) U{x) , for a.e. x G Spec(V) (57) 

where {[B]x}b{x) stands for the matrix representation of [B]x with respect to the ONB B{x) of 
Jv(x); in particular, using that rk(/ — P{x)) = n — d, we conclude that rk([i?]a;) < n — d for a.e. 
X G Spec(V). On the other hand, by construction of B{x), we have that {[-^x^b^x) = D\{x): thus, 
using Eqs (iM|) and (f57)l we have that 

{[A]x + [B]x}b{x) = {[-4]a:}B(x) + {[Bx]}b{x) = U{x)* U{x) , for a.e. x G Spec(V). 

This last fact implies C = A + B ^ U\>{E{Q)) satisfies that Aj([C']a;) = for j G N and a.e. 

X G Spec(V). □ 
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